Pseudogap-related charge dynamics in layered-nickelate i22-xSr^Ni04 (x ~ 1) 
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Charge dynamics and its critical behavior are investigated near the metal-insulator transition of 
layered-nickelate ii2-a;Sra;Ni04 {R = Nd, En). The polarized x-ray absorption spectroscopy exper- 
iment clearly shows the multi-orbital nature which enables the — j/'^-orbital-based checkerboard- 
type charge ordering or correlation to persist up to the critical doping region {x --^ 1). In the 
barely metallic region proximate to the charge- ordered insulating phase, the nominal carrier density 
estimated from the Hall coefficient markedly decreases in accord with development of the pseudo- 
gap structure in the optical conductivity spectrum, while the effective mass is least enhanced. The 
present findings combined with the results of recent angle-resolved photoemission spectroscopy show 
that the pseudogap in the metal-insulator critical state evolves due to the checkerboard-type charge 
correlation to extinguish the coherent-motion carriers in a characteristic momentum (k)-dependent 
manner with lowering temperature. 

PACS numbers: 71.30.-(-h, 78.20.-e, 74.72. Kf, 78. 70. Dm 



I. INTRODUCTION 

Charge ordering is one of the generic phenomena in 
strongly correlated electron systems^ and occurs with a 
variety of ordering patterns (e.g. in forms of stripe and 
checkerboard) depending on the band filling of the con- 
duction electrons. The ordering in the charge sector by 
strong correlation is occasionally accompanied with the 
concomitant or subsequent orders of spin, orbital, and 
lattice degrees of freedom, which enrich the relevant elec- 
tronic phase transitions. One such example is melting of 
the charge order in the critical states, which leads to 
emergence of, for example, the high-temperature super- 
conductivity in layered cuprates^ and the ferromagnetic 
metal in colossal magnetoresistive manganites^. In the 
metallic state derived from the melting of long-range or- 
ders, the charge as well as relevant spin/orbital corre- 
lations may have a large impact on the charge dynam- 
ics, giving rise to rich spatiotcmporal charge responses 
such as formation of pseudogap structure in the momen- 
tum (k)-space. The purpose of this work is to establish 
the link between the spectroscopically-resolved pseudo- 
gap and the charge transport properties in barely metal- 
lic layered nickelates, a typical electron-correlated system 
with charge-ordering instability. 

Layered-nickelate i?2-£cSra;Ni04 {R: rare earth el- 
ements) is a prototypical system showing both the 
charge ordering^ and metal-insulator transition with hole 
dopingi^i^. Besides the diagonal stripe-type commensu- 
rate or incommensurate charge-spin ordering around x = 
1/3^"— the checkerboard (CB)-type charge ordering as 
shown in Fig. 1(a) emerges on the Ni02 plane centered 
around the other fixed point x = \ At a; = 0.5, the 

x^ — y^-orbital-based CB-type charge order develops be- 
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FIG. 1: (Color online). Schematic of electronic struc- 
ture evolution with hole doping in the single-layer nickelate 
i?2-a;Sr3;Ni04. (a) Checkerboard (CB)-type charge ordering 
on the Ni02 plane at a:: = 0.5. (b) Underlying two-dimensional 
hole Fermi surface with k-dependent pseudogap at x = 1.1, 
reproduced from the previous ARPES study.-i^ Solid and dot- 
ted curves represent the Fermi momentum positions where 
one can see a quasiparticle peak and a high-energy pseudogap 
structure, respectively. A three-dimensional electron pocket 
centered only in fc^ = plane is also shown. 



low Tco ~ 500 K under the strong intersite Coulomb re- 
pulsion and cooperative Jahn- Teller distortion. The CB- 
type charge ordered phase has been actually confirmed 
at least up to a; = 0.7, while it has never been observed 
below X = 0.5. In the in-plane optical conductivity spec- 
tra at the ground state, a real-gap structure exists with a 
maximum gap size at x = 0.5,^^ and it gradually evolves 
into a pseudogap in the metallic region proximate to the 
insulating phaseji^ Recently the corresponding structure 
has been found in a k-dependent form by angle-resolved 
photoemission spectroscopy (ARPES) at x = 1.1.— As 
sketched in Fig. 1(b), the main Fermi surface is a (tt, tt)- 
centered large cylindrical hole surface, while the electron- 
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pocket-like Fermi surface with dominant — orbital 
character is formed around the F point. The gap of 
about 0.1 eV is observed to open around the (7r,0) re- 
gion at low temperature. This feature bears analogy 
to the high-energy pseudogap observed in underdoped 
superconducting cuprate o^^i^° and also in colossal mag- 
netoresistive manganites^l. Since the antiferromagnetic 
interaction in the nickelate parent material is one order 
of magnitude smaller than that of the cuprates;^ it has 
been suggested that a plausible origin of the pseudogap 
is not the spin correlation but rather the remnant charge 
correlatiou/ii Here first we discuss the doping variation 
of the hole orbital states based on analysis of linearly po- 
larized x-ray absorption spectra (XAS), and clarify the 
orbital nature of the CB-type charge ordering or correla- 
tion, which can persist up to the metal-insulator critical 
doping region. Then we report on temperature and dop- 
ing variations of the high-energy pseudogap as probed by 
optical spectroscopy and their interrelation with anoma- 
lous charge transport in the course of the metal-insulator 
transition. The present systematic study in the nicke- 
late system can give important hints to understanding 
of generic charge dynamics under the influence of strong 
real-space charge correlation. 



II. EXPERIMENT 

Single-crystalline samples of La2-xSra;Ni04 {x = 0.5), 
Nd2-:,Sr^Ni04 (x = 0.3-1.3), and Eu2-xSr:rNi04 {x = 
0.7-1.2) were grown by a floating-zone method. The 
sample growth above x = 0.8 was carried out in a 
high-pressure oxygen atmosphere up to po^ = 60 atm. 
R = Nd and Eu crystals obtained in the highly-doped 
region were adopted for the following measurements; in 
this system the electronic properties have been reported 
to be least affected by variation of the R speciesj ^°i^^'^^i^^ 

XAS experiments were carried out at undulator beam 
line BL-16A of KEK Photon Factory. The O K-edge 
data were taken in the fluorescence-yield (FY) mode with 
a typical attenuation length of about 600 A. The mea- 
surements were performed in high vacuum of 1 x 10~* 
torr at 50 K. The XAS spectra for E |j a and E || c were 
measured utilizing vertically and horizontally polarized 
soft-x-rays in grazing incidence geometry to the ab plane, 
respectively. The incidence angle relative to the surface 
was set to 10° and the degree of the linear polarization 
was about 95%. Range of the resultant error is shown 
including one caused in the analysis procedure. 

Longitudinal and Hall resistivities and heat capacity 
were measured with the Quantum Design Physical Prop- 
erty Measurement System. In-plane reflectivity spectra 
in the temperature range of 10-600 K were measured 
between 0.01 and 5 eV. The measurements above room 
temperature were performed in a flow of O2 gas to pre- 
vent the reduction, and we conflrmed that the degrada- 
tion at elevated temperature was negligible by thermally- 
repeated measurements. Synchrotron radiation at UV- 



SOR, Institute for Molecular Science, was utilized for 
the measurements (only at room temperature) between 
5 and 40 eV. For Kramers-Kronig analysis, the spectrum 
above 40 eV was extrapolated by a uj~^ function, while 
below 0.01 eV the Hagen-Rubens relation or the constant 
reflectivity was assumed depending on the metallic or in- 
sulating ground state. 



III. RESULTS AND DISCUSSION 

In Fig. 2(a) we show the polarization-dependent O 
K-cdge spectra at a; = 0.5 for E || a and E || c. Aside 
from the main peak (01s — )■ 02p) at about 536 eV, some 
pre-edge peaks assigned to excitations brought about by 
the p — d hybridization (01s — )■ Ni3d)^'2i are clearly 
observed in the lower energy region (3 peaks for E || a and 
2 peaks for E || c). Considering all the possible transition 
processes as sketched in Fig. 2(b), these pre-edge peaks 
at about 532 and 529 eV for x — 0.5 are ascribable to the 
processes A and B, as assigned in the previous worki^ 
The two initial states are Ni^"*" one and Ni^+ one with 
unoccupied x'^ — orbital, respectively, indicating that 
the holes below x = 0.5 are certainly doped into the 
x^ — orbital. Figures 2(c) and 2(d) show the detailed 
doping variation of the O K-edge spectra in the pre-edge 
region for E || a and E || c, respectively. Judging from the 
change below and above a; = 0.5 for E || c, an additional 
peak, which is clearly discerned above x = 0.6, can be 
assigned to the process C, indicating the existence of Ni'^"'" 
site with unoccupied 3z^ — orbital. Above x = 1.0 the 
total peak positions largely shift to lower energy both for 
E II a and E || c; this may arise from the appearance of 
an additional peak, which can be assigned to the process 
D relevant to the high-doping induced Ni^"*" sites. 

By fitting the pre-edge peaks B-D with Gaussian 
curves and comparing the intensity variations for E || c, 
the valence- and orbital-state changes of the Ni site with 
doping X can be roughly estimated. Figure 3(a) summa- 
rizes doping variation of the estimated site occupancy of 
the respective valence and orbital states. Up to x = 0.5, 
the Ni'^+ site with unoccupied x'^ — orbital increases 
nearly linearly in density. This change corresponds to 
the fact that the x"^ — y^-orbital electrons and holes form 
the stripe- and CB-type charge orders on the Ni02 plane. 
Above X = 0.5, on the other hand, the Ni'^"'" site with the 
3z^— hole starts to increase, while that with the x'^ — y'^ 
hole is nearly constant, namely 1/2 per Ni site. This non- 
monotonic variation strongly indicates that the excess 
holes above x = 0.5 are mainly doped into the 3z^ — 
orbital and thus the x^ — y^-orbital-based CB-type charge 
ordering or correlation robustly persist far beyond the 
fixed point a; = 0.5. A similar a;-independent commensu- 
rate order pattern has been reported for a wide doping re- 
gion 0.3 < a; < 0.5 in the so-called CE-type spin-charge- 
orbital ordering of Pri_a;Caa;Mn03. This has been at- 
tributed to occupancy of the other 3z'^ — orbital states 
{z II c— axis) by the extra electrons, which has been con- 
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FIG. 2: (Color online), (a) Polarization-dependent XAS spec- 
tra at O A'-edge for E |1 a and E || c. Some pre-edge peaks as 
labeled with letters (A-D) are clearly discerned in the lower 
energy region (528-533 eV) next to the main peak centered 
at about 536 eV. (b) Schematic of all the transition processes 
corresponding to the pre-edge peaks. Filled arrows added 
in the left and right sides represent the electrons which are 
transferred from the O Is state to the initial Ni 3d eg orbital 
state shown in the middle. Detailed doping variation of the 
O -ft'-pre-edge spectra for (c) E || a and (d) E || c. Gaussian 
fittings for the pre-edge peaks corresponding to the transi- 
tion processes (B, C, and D) as depicted in (b) are indicated 
by dashed, solid, and chain curves, respectively. The peaks 
corresponding to the bracketed final states in the panel C of 
(b) cannot be clearly discerned probably because such higher- 
energy peaks are easily hidden by the main or other pre-edge 
peaks. We confirm that variation of the fitting procedures 
even considering the extra peaks causes negligible changes in 
the estimated site occupancy shown in Fig. 3(a) within the 
error bar and no difference in our conclusion. 



firmed by the electron diffraction measurement Thus, 
the present behavior in the layered nickelate can be re- 
garded as the hole counterpart to the manganite case, as 
also expected in the previous work.^^- The observed ten- 
dency of the hole orbital states is also consistent with the 
recent soft-x-ray- ARPES result suggesting that the holes 
are expected to be nearly equally doped into the two eg 
orbitals at x = 1.1 The existence of nominal Ni''+ site 
above x = l.Q appears to suppress the CB-type ordering 
or correlation and then drive the succeeding insulator- 
metal transition. In addition, in accord with the change 
of the orbital sector, lattice constants as shown in Figs. 
3(b) and 3(c) exhibit the concomitant nonmonotonic de- 
pendence. With sharp increase in density of the Ni'^^ site 
with the 3^^ — r^-orbital hole, the a- and c-axis lattice pa- 
rameters start to increase and decrease above x ^ 0.5, re- 
spectively, reflecting the strong electron-lattice coupling 
characteristic of the eg orbital systems. These exper- 
imental results indeed support that the reported high- 
energy pseudogap is ascribable to the short-range or dy- 
namic CB-type charge correlation persisting due to the 
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FIG. 3: (Color online), (a) Doping variation of the site occu- 
pancy with the respective valence and orbital states for the 
hole doping in Nd2-2,Sr3;Ni04, as estimated by fitting and 
quantifying the peaks B-D in the O K-pre-edge spectra for E 
II c. Room-temperature lattice constants along (b) a- and (c) 
c-axis directions for Nd2-iSriNi04 single crystals. The crys- 
tal structure is low-temperature orthorhombic for x < 0.45 
and high-temperature tetragonal for x > 0.50, respectively. 



orbital degree of freedom. 

Figure 4(a) gives the doping variation of the in-plane 
resistivity pab- Pab gradually decreases with further 
doping holes into the CB-type charge ordered phase. 
Above a; = 1.0 the conductivity (7abi= ^/Pab) shows 
a finite value do in the zero temperature limit, which 
indicates that the metallic ground state emerges. The 
low-temperature weak upturn therein can be fitted with 
the equation aab = ctq + ^-s/T /^'^^ suggesting three- 
dimensional weak localizationj^l The pseudogap forma- 
tion cannot be accounted for in terms of the so-called 
soft-gap,^ - since the energy scale of such weak localiza- 
tion should be far smaller than that of the high-energy 
pseudogap, while that might be an origin of the dis- 
crepancy between the conductivity and the extrapolated 
value in the optical conductivity spectra shown later. 

In Fig. 4(b) we show the temperature and doping 
dependence of the in-plane Hall coefficient i?H- In the 
metallic region at x = 1.0-1.3, positive values nearly 
saturated in the high-temperature limit (^ 4 x 10~^ 
cm^/G) correspond to the case of one hole carrier per 
Ni site, indicating a typical metal with the large Fermi 
surface. Namely, the (tt, 7r)-centered large hole surface 
(Fig. 1(b)) is interpreted as dominantly contributing to 
i?H and its temperature variation, regardless of the ex- 
istence of the electron pocket.— For a: = 1.0 and 1.1, 
i?H shows monotonic increase with lowering tempera- 
ture, while it becomes nearly-temperature independent 
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FIG. 4: (Color online). Temperature and doping variations 
of the in-plane (a) resistivity pat and (b) Hall coefficient Rh 
in J?2-iSra;Ni04. Dotted, solid, and broken lines show the 
results obtained for R = La, Nd, and Eu, respectively. Little 
R dependence is observed as also confirmed in the previous 
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(c) Doping dependence of the low-temperature 



specific heat in a form of C /T vs. T plot for R — Eu. 
Straight lines represent the extrapolation toward zero tem- 
perature to estimate the electronic specific heat coefficient 
(7, Fig. 6(c)). The deviation from the lines below T ~ f 
K arises from the nuclear Schottky component. In addition, 
the curved shape with the steeper slope at a; = 0.7 is due to 
an additional low-temperature entropy component appearing 
around the low-doping region/^ 
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FIG. 5: (Color online), (a)-(c) Temperature variation of the 
in-plane optical conductivity spectra (Tai,(a;) in Nd2-a;Sra;Ni04 
with X = 1.0, 1.1, and 1.3. (d)-(f) Temperature depen- 
dence of the inverse of the Hall coefficient 1/7?h and the low- 
energy spectral weight (effective number of electrons AA'cfi 
at Lo^ = 0.2 eV, see text for details) for the corresponding 
doping levels. Remarkable decrease and increase of 1 /_Rh and 
AA'^off are observed below the characteristic temperature T* 
(indicated by vertical shadows), respectively. 



at a; = 1.3. As discussed later, the pseudogap formation 
is a plausible origin of this anomalous metallic behav- 
ior. Below X = 0.9, on the other hand, the sign of i?H 
becomes negative in the high-temperature region. This 
doping-induced sign change results from the expansion 
(shrink) of the observed electron (hole) Fermi surfaces 
with increasing the electron number. In other words, 
the sign change of i?H in the high-temperature metallic 
state also confirms the multi- band/orbital nature as re- 
vealed by the XAS study. The low-temperature rapid 
divergence with sign inversion is therefore ascribable to 
topology change of the predominant Fermi surface from 
an electron surface to a small hole pocket through the 
antiferromagnetic spin ordering (Tn ^ 80 K).^^ 

Figure 4(c) presents the doping variation of the low- 
temperature specific heat in the layered nickelates with 
no magnetic moment on the i? = Eu site. In the C/T 
vs. plot, the electronic specific-heat coefficient 7 

can be estimated as an intersection of the extrapola- 
tion toward zero temperature. Above x = 1.0, 7 cer- 
tainly has a finite value as a sign of the metallic ground 
state. A noteworthy feature is that its magnitude re- 
mains small (~ 6 mJ/K^ mole) and shows no enhance- 
ment in the vicinity of the metal-insulator transition. 
This behavior is in contrast with the critical mass renor- 



malization effect, which is as intrinsically observed in 
typical filling-control Mott transition systems such as 
Lai.^jSr^rTiOa^o and Lai_:rSr^V0;j2^ The resuh sug- 
gests that the metal-insulator transition in the layered 
nickelate is driven not by the carrier mass divergence 
but by the carrier number vanishment. This situation is 
analogous to the superconducting layered cuprate case as 
studied in La2-a;Sra;Cu04, in which the estimated 7 value 
in the normal state shows no critical enhancement (~ 9 
mJ/K^ mole)'^^'— as a hallmark of the carrier- vanishing 
Mott transition under the presence of the pseudogapiiiiS 
Now let us move on to the detailed charge dynamics 
in the metallic region adjacent to the insulating phase. 
Figures 5(a)-(c) show the in-plane optical conductivity 
spectra ijab{y>) ranging from the anomalous metallic state 
(x — 1.0-1.1) to the normal metal (x = 1.3). We first 
focus on the ground state at 10 K. With increasing x, 
a mid-infrared peak with the real-gap structure (~ 0.6 
eV at a; = gradually evolves into the pseu- 

dogap structure at a; = 1.0-1.1 and then turns to the 
normal metallic spectrum with increasing conductivity 
toward zero energy at x — 1.3. This doping variation 
can be explained in terms of the CB-type charge cor- 
relation persisting far beyond x = 0.5. The pseudogap 
magnitude (~ 0.2 eV) provides a scale of the binding 
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energy of the CB-type-correlated holes, as determined 
by the intersite Coulomb and breathing-type Jahn- Teller 
interactions. As suggested for the case of the layered 
manganites,^^ softening of such a breathing phonon mode 
can potentially drive the pseudogap formation through 
the metal-insulator transition. With increasing x above 
X = 0.5, the gap magnitude decreases reflecting the grad- 
ual weakening of the charge-order amplitude, and the 
remnant polaronic behavior may give rise to the incoher- 
ent carrier dynamics in the anomalous metallic region. 
While little temperature dependence is observed for the 
Drude-like response at x = 1.3, the pseudogap structure 
at X = 1.0-1.1 easily collapses into a broadened gapless 
feature with increasing temperature, and shows satura- 
tion of the temperature dependence at several hundred 
kelvin. Such a spectral variation indicates that the CB- 
type charge correlation develops below some characteris- 
tic temperature in the pseudogap regime. 

To analyze the temperature-dependent pseudogap, we 
show in Figs. 5(d)-(f) the temperature variation of the 
spectral weight below the pseudogap energy, namely the 
effective number of electrons AA^off at 0.2 eV. Nes is 
defined by the following relation. 



2moV 



a{uj)duj. 



Here mo, V , and ojc (= 0.2 eV) are the free-electron mass, 
the cell volume containing one Ni atom, and the cut-off 
energy determined from the isosbestic point energy of 
aab{^)- AAefF at a given temperature T is defined as 
AA"off — Acff(600K) — Aoff(T), representing the amount 
of the spectral weight loss below 0.2 eV. For x = 1.0 
and 1.1, the AN^s value starts to notably increase below 
T* ^ 400 and 250 K, respectively. The characteristic 
temperature T* can be interpreted as an onset tempera- 
ture of the pseudogap formation. Together with AN^ff, 
we show the Hall coefficient in a form of 1/-Rh, which 
also provides a measure for the carrier number in the 
present case of the nearly single band. Below T*, l/i?H 
at X — 1.0-1.1 shows remarkable decrease from the high- 
temperature value corresponding to the large hole surface 
volume without the pseudogap. At x = 1.3, on the other 
hand, the decrease is hardly observed. A plausible origin 
of the temperature variation is the partial disappearance 
of the Fermi surface or the decrease of the mobile carrier 
number due to the CB-type charge correlation. Our re- 
sults indicate that in the pseudogap regime the coherent 
carrier dynamics is partly suppressed in response to the 
strength of the evolving charge correlation. In particular, 
its suppression is expected to selectively occur centered 
around the (7r,0) region (Fig.l (b)), mutually connected 
with the CB-type (tt, tt) wave vector. 

In Fig. 6 we summarize temperature and doping 
variations of the various physical quantities character- 
izing the metal-insulator transition in the layered nick- 
elates. Figures 6(a) and 6(b) give the evolution of dab 
and i?H in the T-x phase diagram, respectively. The 
(Tab downturn temperature systematically decreases along 
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FIG. 6: (Color online). Evolution of (a) the in-plane conduc- 
tivity (Tab and (b) the Hall coefficient 7?h in the T-x diagram. 
Doping variation of (c) the electronic specific-heat coeflicient 
7 and (d) the pseudogap (PG) magnitude Zipc (see text for 
details). They indicate the insulator (I)-metal (M) ground 
state boundary and the proximate critical doping region (CR) , 
roughly denoted by the vertical broken lines, (e) Electronic 
phase diagram of 7?2-a;Sra;Ni04 around the metal-insulator 
transition. The pseudogap formation temperature T* (see 
text for details) are presented together with the checkerboard 
(CB)-type charge order and antiferromagnetic (AF) spin or- 
der transition temperatures (upward triangles;^- downward 
triangles and squares^^). The open and filled marks indi- 
cate the data for i? = La and Nd, respectively; in this system 
the critical temperatures of the spin and charge orderings have 
been confirmed to be least infiuenced by the variation of the 
R species.— 'i^'i^ The solid and broken curves are merely the 
guide to the eyes. 



with Tco, and the metallic ground state emerges above 
X — 1.0 as confirmed by the finite value of a at at the 
lowest temperature. As indicated in the sign inversion 
of the high-temperature i?H value, the dominant Fermi 
surface changes between x = 0.9 and 1.0 from the F- 
centered electron one to the (tt, 7r)-centered hole one. In 
the anomalous metallic state {x = 1.0-1.1), i?H shows 
large increase below T*, while the nearly temperature- 
independent behavior is observed in the normal metal 
(x — 1.3). In addition, 7 in the metallic region is rel- 
atively small and shows no critical enhancement at the 
boundary of the metal-insulator transition, as seen in Fig. 
6(c). These findings again demonstrate that the metal- 
insulator transition in this nickelate system is driven not 
by the carrier mass divergence but by the carrier num- 
ber vanishment due to the pseudogap formation. In Figs. 
6(d) and 6(e) we show the doping variation of the pseudo- 
gap magnitude Apq estimated from the isosbestic point 
of (Tah(a;), together with T* , Tco, and Tn- With increas- 
ing the hole doping level above x = 0.5, Zipc and Tco 
gradually decrease in contrast to Tn , reflecting the weak- 
ening CB-type charge correlation. Even in the critical 
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region at x = 1.0-1.1, however, the pseudogap structure 
with Z\pG of about 0.2 eV develops with lowering tem- 
perature below T*, showing similar temperature-doping 
dependence as in the layered cuprates. Combined with 
the former ARPES study,.^^ it is likely that below T* 
the coherent carrier dynamics around the (tt, 0) region is 
strongly and momentum (k)-selectively suppressed with 
evolution of the CB-type charge correlation. In the nor- 
mal metal at a; = 1.3, the coherent state is expected to 
appear all over the momentum region of the large hole 
surface. 



IV. CONCLUSION 

In conclusion, first we have reported on the detailed 
doping variation of the hole orbital states in the layered 
nickelate, and confirmed the possibility that the CB-type 
charge correlation can persist up to the critical doping re- 
gion because of the multi-orbital nature. Then we have 
systematically investigated the pseudogap dynamics and 
concomitant critical behavior in the metal-insulator tran- 
sition. While no critical mass enhancement is observed 
at the boundary of the metal-insulator transition, the 



carrier number markedly decreases in accord with clear 
evolution of the pseudogap structure. Our results indi- 
cate that the CB-type charge correlation indeed domi- 
nates the charge dynamics in the pseudogap regime and 
that the pseudogap evolution induces the metal-insulator 
transition by gradually and momentum (k)-selectively 
suppressing the density of states or the quasiparticle 
spectral weight at the Fermi level. The present findings 
will provide better understanding of charge dynamics 
near the metal-charge-ordered-insulator transition as well 
as of the k-dependent pseudogap formation in a broader 
class of systems. 
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